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1 Introduction 

Compartmental systems are mathematical systems that are frequently used 
in biology and mathematics. Also a subclass of the class of chemical processes 
can be modeled as compartmental systems. A compartmental system consists 
of several compartments with more or less homogeneous amounts of mate- 
rial. The compartments interact by processes of transport and diffusion. The 
dynamics of a compartmental system is derived from mass balance consider- 
ations. 

The mathematical theory of compartmental systems is of major importance: 

it is the bread-and-butter of analysis for medical researchers, pharmacokineti- 

cists, physiologists, ecologists, economists as well as other researchers 

[Cobelli, C, Romanin Jacur, G. (1976)], [Fagarasan, J.T., DiStefano, J. (1987)], 

[Cobelli, C, Lepschy, A., Romanin Jacur, G., Viaro, U. (1984)], [Jacquez, J. A. (1999)], 

[Nestorov, I. A., Aarons, L. J., Arundel, P. A., Rowland, M. (1998)]. 
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Sometimes it is useful to reduce a model to get a new one with a lower di- 
mension. The technique's name is lumping, i.e. reduction of the number of 
variables by grouping them via a linear or nonlinear function. 

The objective of model reduction methods is to obtain a model that can 
describe the response of the original model accurately and efficiently (cf. 
[Wilkinson, S. J., Benson, N., Kell, D. B. (2008)]). 

Our aim here is to give explicitly possible lumped compartmental systems in 
a few important classes, mainly of symmetric structure such as: mamillary 
models, catenary models and circular models. Some classes can be treated in 
full generality, some only under restrictions on the parameters. 

We also show how to lump systems which are only numerically lumpable. 

The structure of our paper is as follows. In Section 2 the formal definitions 
of reactions, compartmental systems, induced kinetic differential equations 
and that of exact linear lumping are given. Next, our symbolic results are 
presented. Section 4 shows a few examples which had to be treated numer- 
ically. Finally, the results are discussed and further goals are set. We men- 
tion that the present work is a continuation of a few simple statements in 
[Brochot, C, Toth, J., Bois, F. (2005)] on the symbohc lumping of a general 
two compartment model. 



2 Fundamental definitions 



2.1 Reaction mechanism, compartmental system 

A chemical reaction mechanism is a set of elementary reactions. Formally, it 
is a system < M.,7l,a,P >, where 

(1) M and 7^ are sets with M and R elements (M, e N), 7^ = {1, 2, . . . , i?} 
and M = {Xi, X2,..., Xm}, 

(2) a and /3 are matrices with non-negative integers, whose names are stoi- 
chiometric coefficients^ and for which 

(a) for all r E TZ, a{., r) ^ P{-,r), 

(b) if a(., r) = a{.,r ) and /?(., r) = /5(., r ), then r = r , 

(c) for all m e Al there exists r E TZ such that either a{m,r) 7^ or 
/3(m, r) 7^ holds. 
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This mechanism can be represented in the form 

M M 

J2 a{m,r)X^ ^ P{m,r)Xm {r e 7^). (1) 

m=l m=l 

The entities on the two sides of the arrow are the reactant and product com- 
plexes, respectively. 

The number max{X)m=i r),r E TZ} is said to be the order of the reaction; 
thus, first order reactions are obtained if Vr G 7^ J2m=i Oi{m,r) < 1. If in a 
first order reaction it is also true that the length J2m=i Pi^^i ^) of the product 
complexes is also less than or equal to 1, then one has a compartmental system. 
These formal mechanisms are of great practical importance, and are applied 
in many areas as mentioned in the introduction. 

Thus, a compartmental system is a reaction mechanism in which the length 
of all the complexes is not more than one. In this case we only have reaction 
steps of the type Xm — ^p, O, O ^ ij^iV £ Al), where O is the 

empty complex. 

A generalized compartmental system is a reaction in which all the complexes 
contain a single species, and all the species are contained in a single complex, 
i.e. it is a reaction consisting of elementary reactions of three types 

ym^m ^ Vp'^p: ym^m^O, O y^X^, {m^pEM), (2) 

and Xyn is the constituent of a single complex only. 

A generalized compartmental system with no inflow and with some outflow is 
strictly half-open, while it is strictly open if it contains inflows and possibly 
outflows. 

Reaction (1) is said to be mass- conserving if there exist positive numbers 
p(2), . . . , p{N) such that for all elementary reactions 

M M 

a(m, r)p(m) = ^ (3(m, r)p(m) (3) 

m=l m=l 

holds. If the atomic structure of the species are not known, it is not trivial to 

decide whether a reaction is mass-conserving or not [Deak, J., Toth, J., Vizvari, B. (1992)], 

[Schuster, S., Hofer, T. (1991)]. 

A generalized compartmental system is mass- conserving if and only if it is 
closed: the empty complex is not present. 
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2.2 Induced kinetic differential equations 

The usual continuous time, continuous state deterministic model (or, induced 
kinetic differential equation) of reaction (1) describing the time evolution of 
the concentrations Cm is the polynomial differential equation 

R M 

Cm = E(/^(^'0 - a{m,r))kr U cfP'^\ (4) 
where kr denotes the rate coefficient, for all r 

The induced kinetic differential equation of a first order reaction is of the form 

c^Ac + b (5) 

with 

dmp ^0 {m^ p) and &rn ^ (m,p e Ai). (6) 

The induced kinetic differential equation of a compartmental system has an 
additional property 

M 

- amm > J2 ^Prn ("^ ^ M) . (7) 

p^m 

Thus, e.g. there is no compartmental system with the induced kinetic differ- 
ential equation x — x ov with x — —0.5a; + 1/, y = —y + x. 

An easy construction proves that the converse of the above statement is also 
true: a linear differential equation (5) fulfilling the requirements (6) and (7) can 
be considered as the induced kinetic differential equation of a compartmental 
system. 

This statement can be generalized to get our next theorem showing that if 
the right hand side of a kinetic differential equation is the sum of univariate 
monomials and if all the variables have the same exponent in all the rows, then 
- if an additional condition is also met and only then - there exists an inducing 
generalized compartmental system to the system of differential equations. 

Theorem 1 There exists an inducing generalized compartmental system of 
M compartments to the system of differential equations 

M 

C-m ^ ' ^mpiCp)^ ~l~ (8) 
p=l 

(where for all m,p E A4, y"^, y^ G N, y"* 7^ y^, ifm ^ p, a^p, bm ^M.) which is 
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(1) closed, if and only if 6„ = 0, -flmm, a^p, d„, G Eq ; Omm = c?ml/"', 

("^j strictly half-open, if and only if ^ 0,-amm,amp,dm e R^;amm < 

dmU ) 3''^) O'mm ^ , 
f5j strictly open, if and only ifbjji, dmm-) Oimp-i dm ^ -^o"; ^mm — d-mV"^ ■) '^m bjn G 

M+, 

where throughout 

M 

p=\ 

Proof. 

A) The induced kinetic differential equation of (2) is 

P P 

Cm^- y'^iCmf^ E + Z/"' E kmp{Cpf , (3 A;om G M+) 
p p 

Cm — ~ y"^{CmY E ^P™ ^™ E ^mp{CpY + ^mOi (3 A^mO ^ I^^) 
p p 

m e {1,2,...,M}; A;mp eM+; y^' e N;p e {0, 1 . . . , 

Comparing the coefficients we get the only if part of the Theorem. 

B) Given (8) we construct a generalized compartmental system (8) as its 
induced kinetic differential equation: 

yPX/-jd]C y^X^^ fXp^O, o'^^y^Xm, (9) 

(m,p e {l,2,...,M},m^p). 

Reaction (9) induces closed, strictly half-open or strictly open reactions, re- 
spectively. 

2.3 Exact linear lumping 

A special class of lumping is exact linear lumping. 

A system c — foc, with /, c n-vectors can be exactly lumped by an n x n real 
constant matrix Q {n < n), called lumping matrix, if for c = Qc we can find 
an fi-function vector / such that c — foc. 

Not every system is exactly lumpable. A sufficient and necessary condition for 
the existence of exact lumping is Qf{c) = Qf{QQc), where Q denotes any of 
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the generalized inverses of Q, i.e. QQ — In, and In is the fixn identity matrix 
[Li, G., Rabitz, H. (1989)]. 

This condition is equivalent to the requirement that the rows of matrix Q span 
an invariant subspace of f'~^{c) for all c, where f'~^{c) denotes the transpose 
of the Jacobian of / at c. Therefore, in order to determine lumping matrices 
Q we need to determine the fixed /''''(c)-invariant subspaces 
[Gohberg, I., Lancaster, P., Rodman, L. (1986)]. 

In the case of linear differential equation (5), the Jacobian matrix is just A, and 
then f'^{c) — . In this situation, fixed invariant subspaces exist, they are 
spanned by eigenvectors, and they correspond to (constant) eigenvalues. So, a 
linear system is always exactly lumpable and any /''''(c)-invariant subspaces 
will give a lumping matrix. In this case therefore, we have to calculate the 
eigenvectors of . 

We mention here that, if Q is an n x n lumping matrix and P a nonsingular 
matrix of dimension n, then PQ is also a lumping matrix. 

It is not true that a given system can be lumped arbitrarily. For example 

Q kl J ^2 T-) 



cannot lead to a lumped system of the type 5" ^ / — ? P, for 7 := 7i + 

except in the very special case /c2 = /C4, contrary to 

[Conzelmann, H., Saez-Rodriguez, J., Sauter, T., BuUinger, E., AUgower, F., Gilles, E. D.]. 

The next question is whether the lumped system can have an interpretation 
in terms of reactions (or more specially, in terms of compartments), i.e. is the 
lumped system kinetic? To formulate this criterion we use the notion of the 
generahzed inverse matrix [Rao, C. R. (1973)]. 

Farkas [Farkas, Gy. (1999)] gave a sufficient and necessary condition under 
which certain lumping schemes preserve the kinetic structure of the original 
system: A nonnegative lumping matrix leads to a kinetic differential equation 
if and only if it has a nonnegative generalized inverse. 

For the absence of a nonnegative generalized inverse he proved the following 
result: A nonnegative matrix has no nonnegative generalized inverse if and 
only if it has a row such that in the column of each positive entry there exists 
another positive entry. 
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3 Symbolic results 



3. 1 Chains 



In a chain or catenary system the M compartments are arranged in a hnear 
array such that every compartment exchanges material only with its imme- 
diate neighbors and the possible steps are indicated by nonnegative reaction 
rates. The coefficient matrix for a catenary system has nonzero entries only in 
the main diagonal and the first sub-diagonal and in the first super-diagonal. 
The latter case holds only if it is reversible (bidirectional) . 



3.1.1 Irreversible chains 

Let us consider a compartmental system, such as the one in Fig. 2, i.e. a chain 
with unidirectional steps. 




Fig. 1. Irreversible catenary system 



In this case the coefficient matrix A on the right hand side of (5) takes the 
form 

' -k^ ... O" 
A;i -A;2 ... 
k2 -h ... 

. . . kM-2 -ku-i 
... kM-i 

The eigenvalues of the transpose of this triangular matrix are obviously the 
elements on the diagonal: — fci, — /c2, —k^i ■ ■ ■ , —ku-i and (and are the same 
as the eigenvalues of the original matrix). So, the corresponding eigenvectors 



8 



can be found easily, and they take the form: 



[1 
[1 



k\ — k2 
ki 

r-i fci-fca 






(fcl-fc3)(fc2-fc3) 
klk2 









0] 
0] 
0] 



[1 
[1 



fcl-fcM-1 (fcl-fcM-l)(fe2-fcM-l) 



ki 



k\k2 



(fcl— fcM-l)(fc2— fcM-l)---(fcM-2— fcM-l) 
k\k2...kM-l 



0] 

1] 



(Here we only consider the robust case when all the reaction rate coefficients 
are different. Then, the above eigenvectors are independent.) 

If we do not neglect inflows and outflows in a catenary system, the principal 
diagonal of matrix A will change, i.e. instead of —k^ we will have —ki — /li in 
the first M — 1 places (where fii denotes the outflow coefficient for the species 
Xi), and —fiM in the last one, instead of 0. The transpose of the modified 
matrix has the following eigenvectors: 



[1 

fci 

A;i-fc2+//i-/i2 

fclfc2 

(A;i-fe3+/:ii-/:i3)(fe2-fe3+M2-M3) 



[1 






.. 


■ 0] 


1 


.. 


• 0] 


k2 

k2-k3+H2-IJ.3 


1 . . 


• 0] 



1 1 ... 1] , 



corresponding to the eigenvalues —ki — /^i, —k2 — Ai2, ... — kM-i — A*m-i, — A^m- 

The graphical representation in this case, when outflows and inflows are in- 
corporated into an irreversible chain, is: 




Fig. 2. Irreversible catenary system with inflows and outflows 

To get a lumped system for this model, we can take some of the eigenvectors 
above to generate several lumping matrices. 

For example, let us consider an irreversible chain with five compartments. 
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Then the induced kinetic differential equation has the following coefficient 
matrix: 



-ki - Hi 














ki 


-k2 - H2 














k2 


-ks - Hs 














ks 


—ki — Hi 














ki 


-A«5 



Let us compose Q e.g. putting the eigenvectors 

kik2k3 ^'2^3 



fcl-fc2+/il-/i2 



10 



and 



(fel-A:4+/il-/i4)(fc2-A:4+M2-A'4)(fc3-fc4+/i3-M4) ka-kA+fJ-S-fJ-i (fc2-fc4+A»2-M4)(fc3-fc4+M3-M4) 

into it as rows. Then, 



1 



fei-A;2+Mi-M2 (A;i-fc4+Ati-M4)(A;2-A;4+/i2-/i4)(fe3-fe4+/i3-M4) 
fc2A:3 









(A;2-fe4+M2-M4)(fe3-fe4+M3-A'4) 

fe3 

^3-^4+^3-/^4 

1 





After some calculations, we get the lumped system A — QAQ, which induces 
the differential equation below: 



Xi 




— k2 — H2 




Xi 


X2 




—ki — Hi 




X2 



SO we got a new compartmental system with two compartments, where 
ki 



Xi 
X2 = 
+ 



^1 - ^2 + A*i - A*2 



Xi + X2 



kik2k3 



-Xi + 



(ki -k4 + Hi - l^i) (^2 -k4 + H2- fJ'i) {ks -k4 + Hs- l^i) ' 
k2kz , k^ 

■X2 + 1 ; — \ X2, + X4. 



{k2 -k4 + H2- I^4){k3 -k4 + H3- 1^4) ks-k4 + H3- A«4' 



The corresponding reaction (actually, a chain with no interaction between the 
compartments) can be illustrated as follows: 
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or it can be the mamillary system Xi 
induced kinetic differential equation. 



^3 ^2, with ^3 neglected in the 



3.1.2 Irreversible chains with nonuniform directions 

Let us mention here that the irreversible case with nonuniform direction of 
the arrows is simpler than the case of a reversible chain. 

As an example, consider the compartmental system with the following dia- 
gram: 







X 




X 



Fig. 3. An irreversible chain with nonuniform directions 



We can associate to it the kinetic differential equation x — Ax, where 



A 



-ho 

ki k2 

-A;2 - A;3 

k4 

-k4 



The eigenvalues are: —ki, — /c2 — k^, —k^, and 0, with multiplicity 2. From the 

. \k2 + k^ k2 + k^, 
corresponding eigenvectors, [1, 0, 0, 0, OJ, [0, 0, 1, 0, OJ, [0, 0, 0, 0, IJ, 



and 



^3 

-— — - Oil 

k2 k2 



k2 k2 

we can determine a lot of lumping matrices. Depending 



,1,0,0 



on our choice, the lumped system can be kinetic or not. 



For example, if we take 



Q 



1 

10 

1 
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then for the lumped system we get the kinetic differential equation system 



Xi 




-h 




Xi 






-k2 -k3 






X3 




-ki 




X3 



which can be illustrated via the diagram: 



^1 



fc4 



t k2 + k3 



On the other hand, if we take 



Q 



1 0000 
0100 

k2 + k3 k2 + ks 



ko 



ko 



this leads to the matrix 



A^QAQ 



-2ki 
-k2 

2ki{k2 + k3) 

k2 



1 



kik2 



^2 



k2 + k. 



k2 + kj, 
k2 



-ko 



kik2 



k2 + k3 



- k2- k3 



k2+^h 

^2 
k2 + k^ 



ki{k2 + k3) + kl ki{k2 + k3) + k1 



k2 + ks 



It can be seen that in this case the positivity conditions relative to the conve- 
nient elements of the matrix are not fulfilled, as expected in accordance with 
Lemma 1 in [Farkas, Gy. (1999)]. Consequently, A does not result in a lumped 
system which has a kinetic differential equation. The new variables are: 



Xi — Xi 
X2 = Xj, 

k2 + /cs 



2^3 



kr, 



k2 + k3 

Xi H ; X2 + X3. 

ko 



3.1.3 Reversible chains 

To compute the eigenvectors even for a reversible chain consisting of only five 
compartments is unsolvable symbolically. We shall give a numerical example 
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Fig. 4. Reversible Chain 

in section 4.1 below. 



3.2 Mamillary systems 

In these systems all the compartments communicate only with a central com- 
partment, Xm+1-i and there is no direct communication between the other com- 
partments. The possible steps are indicated by nonnegative reaction rates. We 
shall call Xm+\ as the mother compartment and all the other compartments 
will be called daughter or peripheral compartments. 




Fig. 5. Mamillary system 



Only the irreversible case can be treated symbolically; a reversible example will 
be treated numerically in section 4.2. A class of reversible mamillary systems 
with a special structure can still be treated symbolically, this will be shown 
in subsection 3.2.2. 



3.2.1 Irreversible mamillary systems 
Inw£trd flows 

Let us consider an irreversible mamillary system with inward flows such as 
the one in Fig. 6. 
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Fig. 6. A mamillary system with inward flows only 
The coefficient matrix of the reaction rate constants is 

-A;i ... 
-A;2 . . . 

... -kM 

ki k2 ... kM-i kM 

The eigenvalues of the transpose of this lower triangular matrix are obviously 
the elements on the diagonal: —ki, —k2, —k^, . . . , —kM and 0, with the corre- 
sponding eigenvectors: 

[1,0,. ..,0,0], [0,1,. ..,0,0],. ..,[0,0,. ..,1,0], [1,1,. ..,1,1]. 

Notice that if we denote by e, the i-th element of the standard basis for M^, 
i G {1, 2, ... , A^}, then 61,62, ... ,eM create the ffist M eigenvectors of such a 
compartmental system. 

To lump the system of differential equations induced by this model, we can 
choose some of these eigenvectors to generate several lumping matrices. 

In the ffist case, if we do not use the vector [1,1,. . . ,1,1] to generate Q, only M 
of the ffist M elements of the standard basis for M^^"^^, that appear above, we 
will receive a new compartmental system, with M compartments, where the 
new species are taken from the old external ones only. In this case lumping 
actually discards some peripheral compartments and permutes the remaining 
ones. 

If we take an M x M nonsingular matrix, P, i.e. a basis transformation ma- 
trix, then PQ will be another lumping matrix. It will consist of some of P's 
columns, and values being elsewhere. Accordingly the new compartments 
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will be the linear combinations of certain old peripheral compartments. An 
obvious interpretation is that they are measured together. 

Assume e.g. we have chosen Q in the following way: it consists of Cj, ej and 
ek of the natural basis i,j,k e {1,2,...,M}. Let P e R^^^ be an 

invertible matrix. Then 



X = PQx = XiP_i + XjP,2 + XkP.3, 



where p.i-iP.2-,P.3 are the linearly independent columns of P, and the coordi- 
nates of the new composition vector x are linear combinations of the external 
species Xi,Xj and Xk with (in general) different coefficients. 

Now, suppose, the eigenvector [1, 1, . . . , 1] is contained in the rows of matrix Q. 
In this case the system of equation Xi = Y.f=i^ Qij^j defines new compartments, 
composed by some of the existing peripheral ones, plus the sum of the original 
one. 



As an example let us consider the following irreversible mamillary system with 
inward fiows: 



Ik. 

V ^\ V -vj 
— > ^4 < — ^2 



The induced kinetic differential equation is i; = Ax, where 



-h 

-^2 
-ksO 
ki k2 ks 



Using the fact that the eigenvectors of are [1, 0, 0, 0], [0, 1, 0, 0], [0, 0, 1, 0], 
and [1, 1, 1, 1], we can set e.g. 



10 
10 

1111 
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The new variables become 







10 


X2 




10 


X3 




1111 



and the lumped system has the variables 



Xi 
Xi 



Xi = X3 
X2 = Xi 

Xs = Xi+ X2+ Xs + X4. 



The resultant process obeys a differential equation 

-h 
-/ci 




Xi 




X2 




X3 







Xi 




X2 




X3 



which is the induced kinetic differential equation e.g. of the reaction 



^^h3,oh X2, 

that is no more a mamillary system. Or, again, we can take Xi ^ ^ X2, 
and say we are not interested in the change of concentration of X^; we consider 
it as an external species. 

Outward flows 




Fig. 7. Outward flows 
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The induced kinetic differential equation is i; = Ax, wfiere 



... ki 
... k2 



... kM 
...0 -K 



witfi K = (/ci + /c2 + • • • + /cm)- Tfie transpose of it, A^, has a single eigenvalue 
—K with the eigenvector [0, 0, . . . , 0, 1], and an eigenvalue with multiplicity 
M, with the corresponding independent eigenvectors 



[10.. 
[0 1.. 



kM 

kM 



0] 
0] 



[0 0.. 
[0 0.. 



fcM- 



kM 

X 

kM 



0] 

1] 



If we build up a lumping matrix, Q, we get reasonable result only with eigen- 
vectors belonging to the multiple eigenvalue 0, since we get A = in all other 
cases, and it is not worth taking such a Q. 

If the eigenvector [0, 0, ... , 0, 1] appears in the lumping matrix, we obtain a 
lumped system, whose coefficient matrix consists of the elements, except a 
single element on the principal diagonal, which has the value —K. This can 
be represented by the extremely simple reaction X ^ O. 



Irreversible mamillary systems with inward and outward flows 



Instead of giving a general treatment we shall take an example again, as in 
subsection 3.1.2. Let us consider the mamillary system below. 
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The coefficient matrix of tfie induced kinetic differential equation is 



-ki 
-A;2 









^3 



-{ks + k4 + h) 



with the simple eigenvalues —ki.^k2,—K := —{k^ + ^4 + ^5), and with 
the triple eigenvalue 0. The corresponding eigenvectors of its transpose are 



[1,0, 0,0, 0,0], [0,1, 0,0, 0,0], [0,0, 0,0, 0,1], 

0,0,1,-^,0,0 . 
ki 



0,0,1,0,0, 



K 



h 

0,0,l,0,--i,0 
k5 



Taking the lumping matrix 



ki 

1 

K 

ks 

1 --^ 
ki 



10 
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we get the lumped system 



xi^xs + —xe 

h 

X2=X2,- -T-X^ 
K4 

X3 = Xi. 



(12) 



In this case the lumped system's differential equation will be very simple: 



Xi 












Xl 


X2 












X2 


X3 







-h 




X3 



We can associate it to the reaction O. 



3.2.2 SimpUcial compartmental systems 

Suppose we have the following formal reaction steps as follows 



0<^X, "L! XjSo, i0Ti<j;i,je{l,2,...,M}. 

CM-j+i 

The fact that the reaction rate coefficients are the same for many reaction- 
antireaction pairs may come from the application when the compartments 
are physically separated (by a membrane e.g.) parts of the space. In gen- 
eral, such kinds of assumption are made in cases when diffusion is mod- 
eled by mass transport between homogeneous boxes; such models often arise 
[Shapiro, A., Horn, F. (1979a)], [Shapiro, A., Horn, F. (1979b)]. 



The transpose of the coefficient matrix of the induced kinetic differential 
tion is 

Co Ci C2 . . . Cm-1 
CM-1 Co Ci . . . Cm-2 

— 

^ — Cm-2 Cm-1 Co . . . Cm-3 



equa- 



Cl 



C2 C3 ... Co 



with Co := — (ci -|- C2 + • — h cm-i + d). Such a matrix (for which every row 
is a cyclic permutation of the top row) is called a cyclic, or circulant matrix. 
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Its eigenvalues can be calculated easily [Gray, R. M. (2006)]. (Certainly, A is 
a cyclic matrix, as well.) 



M-l 



m 
m=0 

M-l 



A2 — Cm£l" 



m=0 



M-l 

'^M = 5^ C^^M-l 



m=0 



where £fe := e m' , (A; = 0, 1, . . . , M — 1) are the roots of unity. 
The corresponding eigenvectors are 



[1 1 

[1 £1 



^1 



1] 



^M-1 
^1 



[1 £m-1 



_M-1 
'M-l 



Here we meet a new problem which we will not discuss here further: obviously, 
in the applications one needs real lumped systems. 



To be more concrete, let us consider the special case (studying the problem 
of complex numbers in the special case) of Fig. 9. (This system is a special 
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, thus the eigenvalues of 



Co Ci C2 
C2 Co Ci 
Ci C2 Co 



are Ai = Co + Ci + C2, A2 = Co + Ci£i + C2£i = Cq + Ci- 

1 - 



l + iV3 



and A3 = Co + ci£2 + £2^2 = co + ci- 



■+C2- 



2 

1 + /v^ 



f C2- 



i^/3 



For ci = C2 the corresponding eigenvectors are [1, 1, 1], [—1, 0, 1], [—1, 1, 0]. In 
this case we can construct a few lumping matrices which lead to a new, simpler 
system. 

Furthermore, if Ci ^ C2 we obtain the following eigenvectors: [1, 1, 1], 
b-c\- iV3{b + c) {b + 2c)sign(6 - c) + iVSb 



{2b + c)sign(6 — c) — iV^c (26 + c)sign(6 — c) — iy/Sc J ' 
\b-c\+ iV^{b + c) {b + 2c)sign(6 - c) - iy/^b 



(26 + c)sign(6 — c) + iy/Sc (26 + c)sign(6 — c) + iy/Sc 
The effect of lumping on qualitative properties 



One of the major questions connected with lumping is: how are the qualitative 
properties of the lumped and of the original system connected? We investi- 
gated this problem in a more general setting in [Toth, J. et al., (1997)]; here 
we add a new statement: suppose we lump a system of M compartments with 
a coefficient matrix having real eigenvalues into a compartmental system of M 
compartments. Then, none of the concentration versus time curves can have 
more than M — 2 local extrema [Pota, Gy. (1981)]. 



4 A few numerical examples 

4.1 A reversible chain 

Consider a reversible chain formed by five chemical species, let them be Xi, X2, X^i '^a 
and X^. Let the forward and reverse reaction rates be 

ki — l,k2 — 4, ks — 2,k4 — 1 and A;_i = 2, A;_2 = 4, A;_3 = 1, A;_4 = 2, 

respectively, as it can be seen in the following chemical mechanism: 
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1 4 2 1 

^2 A3 A5 
14 5 2 



We can associate it with the following induced kinetic differential equation 
system: 

( 

X\ = -Xl + X2, 
X2= Xl - 5X2 + 4X3, 

i;3 = 4x2 - 6x3 + 5x4, 
X4 = 2X3 ~ 6x4 + 2X5, 
X5 = X4 — 2x5. 



The eigenvectors of the transpose of its coefficient matrix are collected in the 
rows of the matrix below: 



0.2 



-0.2 



-0.2 



-0.689897 0.069693 0.240408 0.449489 1 



0.289897 -2.869693 4.159591 -4.449489 1 



-0.2 



-0.2 



If we take as a lumping matrix 



Q 



0.289897 -2.869693 4.159591 -4.449489 1 



0.2 



-0.2 



-0.2 



after some calculations we receive 



-10.898979 
-2 



and we obtain the lumped model 
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4-2 A reversible mamillary system 



In the following, consider a reversible compartmental system with five com- 
partments. Let A5 be the mother compartment, and Xi, X^, A3, A4, the periph- 
eral ones. Suppose that all of the reaction rates corresponding to the reactions 
from the mother compartment to the peripheral ones have the same value, K. 
Whereas, the reverse reactions also have identical reaction rates, k. 

To this chemical mechanism we can set up the system 



Xi = 


—kxi + Kx^ 


i2 = 


—kX2 + KX5 




—kxs + Kx5 


X4 = 


—kx4 + Kx^ 


X5 = 

< 


k{xi + X2 + X2, + Xi) - AKx 



which describes the time evolution of the concentrations of the species taking 
part in the reaction. Consequently, the coefficient matrix will be 

-A; K 



-A; K 



-A; K 



-A; 



k k k k -AK 



Its transpose has a triple eigenvalue, —k, and two single eigenvalues, and 
and —k — AK. The corresponding eigenvectors are as follows: [—1,0,0,1,0], 

[-1, 0, 1, 0, 0], [-1, 1, 0, 0, 0], [1, 1, 1, 1, 1] and [- A 
Now, we can take several lumping matrices. For example if 



Q 



-1 
k 











A; 



k 



1 

k 



AK AK AK AK 
-11 
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we obtain the lumped system 



X3 



—kxi 

-{AK + k)x2 
-kx3, 



and the corresponding model is 



^ 4K + k 
X2 



4.3 Cycles 



4-3.1 Irreversible cycles 




Fig. 10. Irreversible circular system 

Consider an irreversible circular system with three compartments, A*!, ^"2, A3 
and the corresponding reaction rates ki,k2,k3- Then the coefficient matrix of 
the induced kinetic differential equation is 



A 



-ki ks 
ki -k2 
ko -fe 



Since the eigenvectors of are [1, 1, 1], 
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ki + k2-h + ^kl + {k2 - hf - 2ki{k2 + h) 

2fe ' 



k2[- ki + k2 - k^ + ^kf + {k2 - k^y - 2ki{k2 + k^) 

2kik3 



and 



ki-k2 + k3 + ^kf + {k2 - k^Y - 2ki{k2 + k^ 
2k. 



k2[ki -k2 + k3 + Jkf + {k2 - k^f - 2ki{k2 + ks 



2kik3 

ing up Q from the first two eigenvectors, in the specia 

1 

and /c3 = 3 we obtain the lumping matrix Q 



respectively, build- 



case ki = l,k2 = 2 
1 1 
iV2 -2 + 2iV2 



. After 



some calculations we get A 





-3-V2i 

A, we should take, for example ki = 1, k2 = 1/2 and k^ = 5/128. 



3 3 
. To receive a real valued matrix 



We can also illustrate the region of those values ^2 £ [0,20] and k^ G [0,20], 
for which ki = 1 results in a lumped system with kinetic structure, that is, a 
real valued matrix, A (see fig. 11). 



20 - 



Out[69]= 




Fig. 11. The preferred values of ^2 and ^3 are those outside the curve 
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4-3.2 Reversible cycles 



Consider the reversible cycle with five compartments Xi, X2, X^, X4, and A's. 
Suppose the reaction rates are all equal to a positive real number k. 

Then, we can assign to this mechanism a linear difi^erential equation to describe 
the time evolution of the species's concentrations, with coefficient matrix 



2k 


k 








k 


k 


-2 k 


k 











k 


-2 k 


k 











k 


-2 k 


k 


k 








k 


-2 k 



This is a special circular matrix. = A has two double eigenvalues. 



5 + ^5 



and 



5- \/5 



, and a single one, 0. With the corresponding eigenvectors. 



-V5 + 1 
2 ' 2 
-V5-1 V5 + 1 



,-1,0,1 

,-1,0,1 



-v^+1 



^/5 + l -V5-1 
1,— ^ — , — ,1,0 



, and [1,1,1,1, 1] we 



2 ' 2 

can determine several invariant subspaces in order to find lumping matrices. 
Choose, for example. 



Q 



1 1 



_1 ^V|^ 1 



In this case the lumped system will be 





Xi 








X3 





-^k 







5+V5I 





Xi 




X2 




X3 



We can associate it to the model X2 



5+V5 ■ 5+V5 , ^ 
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5 Discussion, plans 



The most important classes of compartmental systems have been reviewed 
from the point of view of symbolic lumpability. Practically interesting lumped 
systems mainly arise from numerical calculations, which can be carried out 
in all cases without difficulties. Wc used the sentence "which is the induced 
kinetic differential equation of the reaction" recurrently. However, given a 
kinetic differential equation the inducing reaction is by far not unique 
[Erdi, P., Toth, J. (1978), pages 67-69]. 



6 Appendix 



Suppose we are given two natural numbers, n and n, n < n, and an h x n 

matrix Q of full rank with real elements. The question arises: what are the 
necessary and sufficient conditions for the existence of a nonsingular n x n 
matrix P such that all elements of PQ are nonnegative ? 

This question is hard enough to answer in a general case. Here is a result when 
n < 2. One can see, if n = 1, the elements of Q must have identical sign, for 
the existence of such P. 



Now, assume n — 2, and take 



Q 



Furthermore, take 



Then we obtain 



Oil ai2 . . . Oln 

021 0,22 ■ ■ ■ 0,2n 



P21 P22 



(13) 



(14) 



PQ 



OllPll + O21P12 ai2Pll + 022^12 • • • OlnPll + a2nPl2 
ailP21 + a2lP22 ai2P21 + 022^22 ■ ■ ■ ainP21 + a2nP22 



The requirement is that all the elements of the matrix above should be non- 
negative real numbers. This assumption, i.e. the inequalities a^jPn + a2jPi2 > 
0, j = l....,n determine half-planes in the plane (pii,Pi2), passing through 
the origin, i = 1,2. Thus, the problem is to find the cases, when the intersec- 
tions of the corresponding planes (which is in accordance with the first, resp. 
the second row in PQ) are not empty. 
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Examine the columns of the matrix Q. We distinguish 9 cases. In what foUows, 
the symbols +, and — indicate the presence of a positive or a negative number 
in the matrix Q. 




Fig. 12. A graphical representation for the case 3. 



The slope of the line with equation anpu + a2iPi2 = depend on the number 
—aii/a2i. Here this fraction is a positive number. The shaded region in the 
figure represents the region that is excluded from the solution. 

Thus, after the geometrical consideration, we can conclude: if the matrix Q 
contains columns of form [l] and Zl simultaneously, V a, 6 G M, with a and 
b different from at the same time, there does not exist nonsingular, 2x2 
matrix P, which satisfies the requirement PQ > 0. This is the case when 
matrix Q has a pair of columns of type 2 and 6, or 4 and 8, or 3 and 7, or 5 
and 9. In the last two cases the elements could only differ in sign. Henceforth, 
for a shortest notation we will use 26,48,37,59 to point to pair of cases when 
the matrix P does not exist. 
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Finding all the cases when the matrix Q contains three columns which precisely 
exclude together the existence of P, we lean on the geometrical representation 
again. Assume Q does not contain a pair of columns fitting the case described 
above. If we check the three half-plane cases, we get the following result: 247, 
257, 258, 358, 368, 369, 469, 479 and 569 are the cases that exclude each other, 
i.e. the existence of P, by all means. Furthermore, there are other instances for 
the nonexistence of such a P, 259, 347, 357, 359, 367, 378, 379 and 459, but 
in these cases we still have to verify another condition regarding the slopes. 
This fact will be illustrated later in an example. 

In the case of the intersection of four half-planes, assumed that we did not find 
in M columns corresponding to either cases given earlier, we get only one case 
for empty intersection, specifically for 2358, i.e. for a matrix that contains the 
columns [°], [t], [^], [q] together, disregarding the order. 

The cases presented previously exhaust all the cases, when to a 2 x n matrix 
given in (13) we cannot find a 2 x 2 nonsingular matrix, P, so that all the 
elements of PQ are nonnegative. 

Let us take an example to illustrate the problem. Consider the matrix 



Q 



5 2 2 -3 
-2 1-1 



and P as in (14). Then 



pg = 



5pii - 2pi2 2pii 2pii + pi2 -3pii - pi2 

5p21 - 2p22 2p21 2p21 + P22 "3^21 " P22 



whose elements must satisfy the system of inequalities: 



5pii - 2pi2 > 0, 



2pii 

2pii Pl2 
-3pii - pi2 > 0, 
5p2i - 2p22 > 0, 

2P21 

2p21 + P22 



>o, 
>o. 



>o, 
>o, 



(15) 



-3p21 - P22 > 0. 



Consider the first four of them, and give a geometrical representation as in 
figure 13. 
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Fig. 13. Example for nonexistence of P Fig. 14. Example for existence of P 

The picture was created with Maple program, and one can see, that in this 
case the inequality system (15) does not have any solution. This was only to 
be expected, because the matrix contains columns of type 3, 4, 5 and 9, and 
359 is a critical case, since the slope of the line corresponding to case 9 is 
smaller then the one corresponding to 5, i.e. — (— 3)/(— 1) < —2/1. 

Notice that if we take ai3 = 18 in Q, the situation will change, i.e. the exis- 
tence of P will be insured, because in this case the direction of the inequality 
regarding the slopes will change, as it can be seen in Fig. 14. 

Choose a point from the region that indicates the solution (the dark region in 
Fig. 14), such as {pn,pi2) = (1, —5). To determine p2i andp22 one must choose 
their values from the same region, except the case when the point (^21,^22) can 
be found on the line defined by the origin and the point {pn,Pi2) = (1, —5). 
The justification of this statement is as follows. The equation of the line that 
passes through the points (0,0) and (pu^Pu) is y = {pi2/pii)x. If (^21,^22) 
is a point on this line, then it must satisfy the relation puP22 — P12P21 = 0, 
which is equivalent to the condition det P = 0, but P cannot be singular. For 
example we can choose (^21,^22) = (1/2, —4). In this case det P = —3/2. Now, 
check the nonnegativity of the elements of the matrix PQ: 



1 


-5 




5 2 18 -3 




15 2 13 2 


1 

. 2 


-4 




-2 1-1 




21 1 c; 5 
_ 2 ^ ^2. 



Suppose we are given two natural numbers, n and h, h < n, and an h x n 
matrix Q of full rank with complex elements. The question arises, what are 
the necessary and sufficient conditions for the existence of a nonsingular nxn 
matrix P such that all elements of PQ are real? 

Observe that for n = 1, Q must have a special form to find a suitable P to it, 
i.e. we must have either Q = [irrii im2 ■ ■ ■ inin] , with nij G M, j = 1, . . . , n. 
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or all the elements of Q have to be real, but this is a trivial case. 

Consider ft > 1 and Q an n x n matrix, with Q = Qi + iQ2 where Qi, Q2 G 
A4nn(^)- A sufficient condition for the existence of a ri x n matrix P such that 
all elements of PQ are real is: 

(1) Qi and Q2 are nonsingular; 

(2) QJQ2 = QlQi. 

In this case we can choose P — Qj — iQ^. Realize that if the conditions above 
are satisfied, then 

PQ = {Ql - iQl){Qi + iQ2) = QlQi + iQlQ2 - iQlQi + QIQ2 
= <5? + Q2 + KQIQ2 - QlQi) e Mnn{^). 



In what follows suppose that the conditions above are not satisfied and take 

a2n + ihn 



Pih + iqih 
P2h + n2h 

The requirement that all the elements of PQ are real is equivalent to a linear, 
homogeneous system 

h 

^{Pijbjk + Qijajk) =0, A; = l~n, l^l,n, (16) 
i=i 

where the number of unknowns are 2n^ and the number of equations are nn. 

Notice that this system has an interesting property: it can be divided into 
independent subsystems regarding the unknowns. In this case each of these 
subsystems can be solved separately. Moreover, all of the results can be writ- 
ten in identical form, since they have the same coefficient-matrix. Therefore, 



Oil -I- ibu ai2 + ihi2 ■ ■ ■ 

021 + ih2\ 022 + «&22 ■ ■ ■ 
Onl + ibni an2 + ibh2 ■ ■ ■ 

Pii + iqu P12 + iqi2 ■ ■ ■ 
P21 + iq2\ P22 + iq22 ■ ■ ■ 

Phi + iqhi Ph2 + iqh2 ■ ■ ■ 
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taking I — 1, consider and treat only the subsystem with n equations 

h 

^{Pijbjk + qijajk) = 0, A; = l~ri, (17) 

j=i 

with unknowns pii,pi2, ■ ■ ■ ,Pi,n, and qn, qi2, ■ ■ ■ , qi,n- Depending on n and n 
the number of unknowns can be smaller or bigger than the number of equa- 
tions. The coefficient-matrix is 

bii 621 • • • bfii dii 0,21 ■ ■ ■ CLfii 

bi2 &22 ■ ■ ■ bn2 ai2 022 ■ ■ ■ Cin2 



bin ^2n ■ ■ ■ bfin Qi^ a2n ■ ■ ■ O'hn 



(18) 



Denote by d the main determinant for (18). Then, if d 7^ 0, system (17) has 
identical zero solution: (0, . . . , 0) . In this case det(P) = 0, so, P does not 

^ V ' 

2n 

satisfy the nonsingularity requirement. 

For an adequate P we must have = 0. If it holds, we must specify the rank of 
(18). Let us denote it by r. If r = 2n, i.e. the system (17) is determined, and has 
a unique solution. But in this case all the other subsystems mentioned above 
for (16) has the same (constant) solution. Thus, the matrix P will consist of 
identical rows, so we have det(P) = again. 

For r = 2n — 1 the system is indefinite, and the solution has the form 

(q;iCi,Q;2C2, ...,Q;2nC2n), 

where are parameters and q represent constants, i = l,...,2n. Conse- 
quently, P can be given in such a way that its rows will be fc-times for the 
others. Thus, we discover a nonsingular P once more. 

Finally, if r < 2n — 1, then there exists an adequate matrix P, because in this 
case the solution of (17) can be expressed with at least two parameters, and 
this gives possibility for choosing linearly independent lines to P. 

Now, look at an example. Consider 

Q = , and take P 



l + i2 + i4 + 2i2 + 2i 
-1 2i U -2 



Pii + iqii P12 + iqi2 
P21 + iq2\ P22 + iq22 
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In this case the coefficient-matrix (18) is 



10 1-1 
12 2 
2 4 4 
2 2 -2 



and for it r = 2 < 2n — 1. Thus, wc can compute a nonsingular P G A^22(C). 
Since the homogeneous, hncar system with coefficient-matrix above has solu- 
tion {—a + P, —^{a + P), a, P), we can take pn = 0,pi2 = 1, Qn = —1, qi2 = — 1 
and similarly p2i = 2,p22 = 2, Q21 = —3, ^22 = —1- After some calculations we 
obtain 

""o 3 6 



PQ 



3 9 18 8 
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